Variational Methods for Biomolecular Modeling

Structure, function and dynamics of many biomolecular systems can be characterized by the energetic variational principle and the corresponding systems of partial differential equations (PDEs). This principle allows us to focus on the identification of essential energetic components, the optimal parametrization of energies, and the efficient computational implementation of energy variation or minimization. Given the fact that complex biomolecular systems are structurally non-uniform and their interactions occur through contact interfaces, their free energies are associated with various interfaces as well, such as solute-solvent interface, molecular binding interface, lipid domain interface, and membrane surfaces. This fact motivates the inclusion of interface geometry, particular its curvatures, to the parametrization of free energies. Applications of such interface geometry based energetic variational principles are illustrated through three concrete topics: the multiscale modeling of biomolecular electrostatics and solvation that includes the curvature energy of the molecular surface, the formation of microdomains on lipid membrane due to the geometric and molecular mechanics at the lipid interface, and the mean curvature driven protein localization on membrane surfaces. By further implicitly representing the interface using a phase field function over the entire domain, one can simulate the dynamics of the interface and the corresponding energy variation by evolving the phase field function, achieving significant reduction of the number of degrees of freedom and computational complexity. Strategies for improving the efficiency of computational implementations and for extending applications to coarse-graining or multiscale molecular simulations are outlined.Comment: 36 page

Vesicle-Like Biomechanics Governs Important Aspects of Nuclear Geometry in Fission Yeast

It has long been known that during the closed mitosis of many unicellular eukaryotes, including the fission yeast (Schizosaccharomyces pombe), the nuclear envelope remains intact while the nucleus undergoes a remarkable sequence of shape transformations driven by elongation of an intranuclear mitotic spindle whose ends are capped by spindle pole bodies embedded in the nuclear envelope. However, the mechanical basis of these normal cell cycle transformations, and abnormal nuclear shapes caused by intranuclear elongation of microtubules lacking spindle pole bodies, remain unknown. Although there are models describing the shapes of lipid vesicles deformed by elongation of microtubule bundles, there are no models describing normal or abnormal shape changes in the nucleus. We describe here a novel biophysical model of interphase nuclear geometry in fission yeast that accounts for critical aspects of the mechanics of the fission yeast nucleus, including the biophysical properties of lipid bilayers, forces exerted on the nuclear envelope by elongating microtubules, and access to a lipid reservoir, essential for the large increase in nuclear surface area during the cell cycle. We present experimental confirmation of the novel and non-trivial geometries predicted by our model, which has no free parameters. We also use the model to provide insight into the mechanical basis of previously described defects in nuclear division, including abnormal nuclear shapes and loss of nuclear envelope integrity. The model predicts that (i) despite differences in structure and composition, fission yeast nuclei and vesicles with fluid lipid bilayers have common mechanical properties; (ii) the S. pombe nucleus is not lined with any structure with shear resistance, comparable to the nuclear lamina of higher eukaryotes. We validate the model and its predictions by analyzing wild type cells in which ned1 gene overexpression causes elongation of an intranuclear microtubule bundle that deforms the nucleus of interphase cells

Remodelling of the angular collagen fiber distribution in cardiovascular tissues

Understanding collagen fiber remodelling is desired to optimize the mechanical conditioning protocols in tissue-engineering of load-bearing cardiovascular structures. Mathematical models offer strong possibilities to gain insight into the mechanisms and mechanical stimuli involved in these remodelling processes. In this study, a framework is proposed to investigate remodelling of angular collagen fiber distribution in cardiovascular tissues. A structurally based model for collagenous cardiovascular tissues is extended with remodelling laws for the collagen architecture, and the model is subsequently applied to the arterial wall and aortic valve. For the arterial wall, the model predicts the presence of two helically arranged families of collagen fibers. A branching, diverging hammock-type fiber architecture is predicted for the aortic valve. It is expected that the proposed model may be of great potential for the design of improved tissue engineering protocols and may give further insight into the pathophysiology of cardiovascular diseases

Modeling morphological instabilities in lipid membranes with anchored amphiphilic polymers

Anchoring molecules, like amphiphilic polymers, are able to dynamically regulate membrane morphology. Such molecules insert their hydrophobic groups into the bilayer, generating a local membrane curvature. In order to minimize the elastic energy penalty, a dynamic shape instability may occur, as in the case of the curvature-driven pearling instability or the polymer-induced tubulation of lipid vesicles. We review recent works on modeling of such instabilities by means of a mesoscopic dynamic model of the phase-field kind, which take into account the bending energy of lipid bilayers